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Abstract. - This paper proposes a simple setup for introducing an artificial magnetic field for 
neutral atoms in 2D optical lattices. This setup is based on the phenomenon of photon-assisted 
tunneling and involves a low-frequency periodic driving of the optical lattice. This low-frequency 
driving does not affect the electronic structure of the atom and can be easily realized by the same 
means which employed to create the lattice. We also address the problem of detecting this effective 
magnetic field. In particular, we study the center of mass wave-packet dynamics, which is shown 
to exhibit certain features of cyclotron dynamics of a classical charged particle. 



Introduction. — A major motivation of the current research with cold atoms in optical 
lattices is the prospect of simulating solid state physics. However, to have a full access to 
solid state-like phenomena in these systems (which can be considered as artificial crystals) 
artificial electric and magnetic fields must be introduced. In present days experiments an 
external electric field is routinely mimicked by accelerating the optical lattice [IH3] , by using 
the gravitational force 4-6 j and by a combination of gravitational and levitational forces 
[7] . The case of artificial magnetic field is more difficult to achieve in laboratory realizations 
because it requires a setup where the atomic wave function acquires a finite phase when the 
atom tunnels along a closed path on the lattice. 

Such a setup was suggested in the seminal paper [S], where the authors used two inde- 
pendent 2D optical lattices for two different internal atomic states, which are coupled by 
additional Raman lasers. The Rabi transition between internal states induces hopping of 
the atom between nearest lattice sites, where the required phase accumulation is achieved 
by using a special geometry for the Raman beams. This idea was further developed in 
Ref. [9,10 . It was shown that one can also introduce non-Abelian gauge potentials by using 
the Raman-laser technique. For atoms in a harmonic trap (no lattice) the Raman scheme 
with 3 spin states of the F = 1 electronic ground state of 87 Rb atom has been recently 
realized [11]. In this setup the magnitude of the effective magnetic field is defined by the 
gradient of the real magnetic field which splits the F = 1 level into Zeeman sublevels. We 
note that the cited experiment also reported a drawback of the Raman-laser based tech- 
niques. Namely, the spontaneous emission from Raman beams kicks atoms out of the trap, 
thus causing the population decay. 

A different method of introducing artificial magnetic field was proposed in Ref. . Un- 
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like the setups of Ref. [8HTT] it does not rely on the internal atomic structure - the magnetic 
field being mimicked by an oscillating quadrupole potential at frequency u>, together with a 
periodic modulation of the hopping matrix elements. Here, the magnitude of the effective 
field is defined by the strength of the quadrupole potential, measured in units of fouj. Unfor- 
tunately, because of a number of approximations involved, this setup does not ensure exact 
equivalence with a real magnetic field. 

In the present work we suggest another, alternative method of creating an artificial mag- 
netic field for neutral atoms in a 2D optical lattice. The method is based on the phenomenon 
of photon-assisted tunneling, which is well studied in the case of (quasi) ID lattices [3j0[7]. 
It has the great advantage over the above discussed setups that it does not involve Raman 
lasers and, yet, ensures complete equivalence with a charged particle in the real magnetic 
field. 

The problem of creating artificial magnetic fields for neutral atoms is closely related to 
the detection problem. Here one addresses magnetic field effects, which can be observed by 
using detection techniques of cold atoms physics. In particular, it was argued in the already 
cited paper [5] that the effective magnetic field produces a specific interference pattern for 
the atomic density (see Fig. Q] below), which carries information about the field magnitude. 
We revisit here the problem of detecting an artificial magnetic field. In particular, we discuss 
a variety of interference patterns developed by the atomic wave function in the course of 
time and study the center of mass wave-packet dynamics, which exhibits certain features of 
the cyclotron dynamics of a classical particle. 

Photon-assisted tunneling and effective magnetic field. The phenomenon of 
photon assisted tunneling refers to a quantum particle (an atom) in a ID lattice, which is 
subject to DC and AC fields Using the tight-binding approximation (which will be our 
theoretical framework from now on) the system Hamiltonian reads 

H = -- J2(\m + l)(m| + /i.e.) + a[F + F u cos(wt + <j>)] ^ | m)m(m\ , (1) 

m m 

where \m) are the Wannier states, a the lattice period and J the hopping matrix ele- 
ment. For vanishing DC and AC fields the egenfunctions of (TTJ) are extended Bloch waves 
with dispersion relation E(k) = — J cos(cik). If F ^ the eigenfunctions become localized 
Wannier-Stark states with an equidistant spectrum with level spacing AE = dF = hcus- A 
periodic driving of the system at the frequency u which matches the Bloch frequency ujb 
couples these Wannier-Stark states into extended quasienergy states with dispersion rela- 
tion E(k) = — Jcos(aK), where J — JJi{F UJ /F) and J\{x) is the Bessel function of the first 
kind| Moreover, if 

aF > J (2) 

the Bloch period is small in comparison with the characteristic tunneling time and the 
dynamics of the system {T]) can be described in terms of the effective Hamiltonian 

# = -~^(|TO+l)He i0 + /i-c) , J=JMF U /F), (3) 

m 

as already shown in detail in Ref. [14.. The additional parameter 4> in ([3]), which was actually 
overlooked in early studies of photon-assisted tunneling, is the phase difference between the 
Bloch and field oscillations. We mention that a possibility of experimental control over the 
parameter tf> has been recently demonstrated in Ref. [TJ. 

1 See Ref. [3][6j[7] for experimental studies and Ref. I131I14| for relevant theoretical analysis, which also 
includes the case of interacting atoms. 

2 In general case the photon-assisted tunneling takes place at integer ration uj/u = n, where J = 
J 3n(a,Fu /hu). It may also happen at rational ug/w, if the Hamiltonian JT} includes hopping to the next 
to nearest-neighbor sites. 
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Fig. 1: Dynamic of the driven system Q, left, as compared to dynamic of the effective system J5J, 
right. Parameters are J x = 0.4, J y = 2, F = 1, F u = 0.4085F (hence J y = 0.4), and a = 1/6. 

Now we consider the quantum particle in a 2D lattice. The lattice is tilted and driven 
in the y direction (index m) and we assume the phase of the AC field to vary linearly along 
the x direction (index I), 

l,m l,m 

+a^[F + cos(u;t - 2iral)]\l, m)m(U m\ . (4) 

Lrn 

Repeating the above arguments for the ID lattice it is easy to show that this setup realizes 
the effective Hamiltonian 

H=-Y^2(\ l + 1 , m )(^ m \ + h.c)- ^(|I,m+l)He iW + /i.c) , (5) 

/ m 

which coincides with the tight-binding Hamiltonian of a charged particle subject to a mag- 
netic fieldH 

To check validity of the approximation involved in the transition from ([?]) to ([5]) we 
simulate the dynamics of both systems for a finite lattice size, 1 < l,m < L, periodic 
boundary conditions, and uniform initial wave function, \^(t = 0)) = ^i,m\l, m), 

^;,m = i/L. It is easy to show that in this case the wave function remains uniform along 
the y direction. Thus, following Ref. [5], we consider the quantity ni(t) = J2 m \^i,m{t)\ 2 ■ 
Time evolution of the density ni(t) calculated on the basis of Eq. §5§ for J x = J y = 0.4 
and a — 1/6, is shown in the right panel in Fig. Q] It is seen that the density develops a 
periodic pattern with the spatial period given by 1 /a and characteristic time period defined 
by the cyclotron frequency ui c = 2iraJ/h. The left panel in Fig. [T] shows the dynamics of 
the system (j4|) for the same a = 1/6 and J x — 0.4, and the other parameters adjusted to 

3 For a charged partical the parameter a = eBa? /he, where e is the charge and B the magnetic field 
magnitude. 



p-3 



A.R.Kolovsky 



Fig. 2: Scheme of the proposed setup, which involves two crossing running-wave beams. Location 
of the 2D optical lattice created by two standing-wave laser beams (not shown) is indicated by the 
circle. In this setup the magnitude of artificial magnetic field (the parameter a) is defined by the 
angle between the running-wave beams. 



have J y = J x , namely, F = 1, F u = 0.4085^, and J y — 2. A comparison of these results 
indicates that the suggested setup indeed introduces an effective magnetic field with the 
magnitude defined by the phase gradient of the driving force. We note that to observe a 
difference between the figures we intentionally choose a relatively small F, - for a larger F 
the figures would be indistinguishable. 

The question of how one can realize the required driving in a laboratory experiment 
is in turn. A possible scheme, which involves two crossing running-wave beams with the 
wave-vector k, is depicted in Fig. [2] In this scheme we assume that (i) the laser frequencies 
are slightly mismatched by Af2, (ii) the center of the 2D optical lattice is shifted from the 
crossing point by the distance yo, and (iii) there is an additional homogeneous static field 
of magnitude Fq . Then, locally, the optical potential created by two running- wave beams is 
given by the following expression: 



dV(r) 

H int =F (y- y ) + —^(v - Vo) sin 2 (qx - Aflt) , 



(6) 



where V(r) is proportional to the envelope function for the light intensity and q = fcsin#. 
In the tight-binding approximation the potential (J6)) takes the form of the driving term 
in the Hamiltonian (gj) with F = F + V/2, F u = V'/2, u = 2AO, and the parameter 
a = (ka/2n) sin (9. 

Different gauges and phase imprinting. — It is worth stressing that the proposed 
setup (as well as all previously suggested setups) realizes not so much a magnetic field but the 
tight-binding Hamiltonian of a charged particle corresponding to a given vector potential. 
In particular, the Hamiltonian 



H = -— {\l + 1, m){l, m\ + h.c.) 
i 



l)(m|e 



i2nal 



h.c. 



(7) 
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which coincides with ([5]), refers to the Landau gauge A = A(0, x, 0). If we choose the gauge 
A = A(—y, 0, 0), the Hamiltonian reads 

JI = - (|( + 1, m)((, m |e- & ™ + k.) - ^ (|i,m + + ft.c.) , (8) 

i m 

while for the symmetric gauge, A = A(—y/2, x/2, 0), one has 

H = + l,m)(l,m\e- ivam + h.c.) - (|Z,m+ l)(m|e w + h.c.) . (9) 

/ m 

Although these Hamiltonians have the same spectrum, their eigenfuctions are different. Thus 
the Hamiltonians ([7][9]) generate different dynamics. In the other words, when evolving in 
time a given initial state, one gets different interference patterns. These patterns are shown 
in Fig. (3][a-c), where we additionally assume the presence of a harmonic confinement, i.e., 
the Hamiltonians are augmented by the term 

ff7 = |l> a + n* a )|J>m><l,m|. (10) 

As the initial state we choose the ground state of the system for a = 0, which for 7 « 1 
is well approximated by the two-dimensional Gaussian G(l,m) of the width a = (J/7) 1 / 4 . 
In the course of time the interference patterns change with the characteristic time period 
defined by the cyclotron frequency, recovering quasi-periodically the initial state. 

It is interesting to note that the interference patterns (b-c) can be reproduced by the 
system ([7} alone, if the initial state has the imprinted phaseQ 

ip l , m (t = 0) = e- i P lm G(l,m), (II) 

where /3 = 2ira and ft — ira, respectively. Since phase imprinting introduces an interference 
pattern by itself, Fig.[3jd) shows the case a — and (3 — n/6 for the sake of comparison. At 
the first glance the depicted pattern resembles that for a = 1/6. However, more thorough 
inspection of this pattern reveals the square symmetry of the lattice, which becomes more 
transparent for larger times. Beside this, for a — 0, the overall width of the packet oscillates 
in time. This breathing of the packet width is absent if a 7^ 0. It should be also mentioned 
that in this section we focus on the short time dynamics for t <C 2-7r/f2 7 , where i7 7 is 
defined below in Eq. (|13|) . For these times the only effect of the harmonic confinement is a 
modification of the boundary conditions. 

We come back to the original system (|4|). We have checked that the system (j4|) well 
reproduces the interference patterns depicted in Fig. [3ja-c) , if condition ([2]) is satisfied. If 
this condition is violated, the continuous dynamics of the driven system ^ may not coincide 
with dynamics of the target system ([5]). However, even in this case the discrete dynamics 
of CP and (0) do coincide]! 

Cyclotron dynamics. — Observation of the interference patterns discussed in the 
previous section requires a 3> l/a. In this section we analyze the opposite situation, where 
the magnetic length l/a is larger than the wave packet width a. In this case the wave 
function does not develop an interference pattern. Nevertheless, one can detect the artificial 
magnetic field by watching the evolution of the wave-packet center of gravity, which is 
expected to follow the classical trajectory of the system. Needless to say that the classical 

4 In a laboratory experiment phase imprinting Ijllll can be done by exposing the system to a quadrupole 
potential for a given time. 

5 By discrete dynamics we mean the dynamics in terms of the system evolution operator over the driving 
period. 
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Fig. 3: Panels (a-c): Interference patterns at t = 107r for J x = J y = 0.2, a = 1/6, and 7 = 2- 10~ 4 , 
calculated on the basis of the Hamiltonians ®> and ([9]), respectively. Alternatively the figures 
show solution of the Schrodinger equation with the Hamiltonian (O for the initial state with 
/3 — 0, /3 = 27ra, and /3 = net. Panel (d): Interference pattern for a = and /3 — n/6. 



trajectories (as well as the wave-packet dynamics associated with these trajectories) do not 
depend on the gauge]! 

The classical counterpart of the quantum tight-binding Hamiltonian ([7]) reads [17] 

H c i = -J x cosp x - Jy cos(p y + 2nax) + ^(x 2 + y 2 ) , (12) 

where the last term takes into account the parabolic confinement and we set the lattice 
period a to unity. As initial conditions we consider p x = p y — 0, and (x, y) = (xo,yo)- 
In a laboratory experiment one realizes these initial conditions by suddenly shifting the 
center of the lattice by the distance r = (x ,2/o) [15] - For a = this shift induces dipole 
oscillations, where the wave-packet center of gravity oscillates around the lattice origin with 
the frequency fid = y/lhj/h, [15jll6]. Unlike these dipole oscillations, for a ^ the classical 
trajectory encircles the lattice origin (see Fig. HJ. Note that for a ^ the Hamiltonian 
(|12l) is not separable and, hence, trajectories may be regular or chaotic. Regular trajectories 
correspond to the stability islands around the points, where arguments of the cosine functions 
in (fT2|) are multiple of 27r and where one may use the effective mass approximation. These 
trajectories encircle the center of the parabolic lattice in the same direction (defined by the 
sign of a) with the frequency 

fl 7 = j/2irha . (13) 

6 For plane lattices (no harmonic confinement) the wave-packet dynamics of the system has been studied 
in the recent paper [17] . It is also shown there that the semiclassical regime of the wave-packet dynamics 
corresponds to the limit a — ¥ 0, where the cyclotron radius essentially exceeds the lattice period. 
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Fig. 4: Left: 4 different classical trajectories for a = 1/6, J x = Jy = 1 and 7 = 0.01. Compu- 
tational time is t = 200-7T, which approximately coincides with T 7 = 209.447T. Right: Occupation 
probabilities at t = 2007T. Initial wave packet corresponds to the ground state of the system for 
a = 0, shifted by 21 lattice periods to the right. 

Chaotic trajectories have no well defined encircling direction and change it in a random 
way. This observation helps us to understand the result of numerical simulation of the 
wave-packet dynamic shown in the right panel in Fig. 01 Here the main part of the wave 
packet rotates counterclockwise with the frequency (fT5]l . while a fraction of the packet goes 
in the opposite direction. Concluding this section we note that the frequency (I13[) can be 
also obtained quantum-mechanically, - it corresponds to the splitting AE = ?if2 7 of the 
low-energy Landau levels of the systems ([7]|9]) caused by the parabolic term (fl"0|). 

Conclusion. — We have proposed a simple setup for introducing an artificial magnetic 
field for neutral atoms in 2D optical lattices. This setup involves a low-frequency driving of 
the optical lattice, which can be easily realized by the same means which employed to create 
the lattice. Note that such low-frequency driving does not affect the electronic structure of 
the atom. This constitutes the main difference of our setup from the previously suggested 
setups, which involve Raman transitions between atomic spin states to mimic the magnetic 

field [EHIMIBH2Q]. 

In the second part of the paper we have discussed two dynamic-based methods of detect- 
ing the (artificial) magnetic field. Importantly, we take into account harmonic confinement 
which is inevitably present in a laboratory experiment. 

The first method is the observation of interference patterns developed by the atomic wave 
function in the course of time. This method assumes a highly coherent initial state (a Bosc- 
Einstein condensate) with the coherence length exceeding the magnetic period. A variety 
of interference patterns signaling the presence of the effective magnetic field is identified. 

The second method does not require a coherent initial state. It relies on the underlying 
classical dynamics of an atom in a parabolic lattice. For vanishing magnetic field these 
dynamics consist of dipole oscillations with the frequency defined by the strength of the 
parabolic confinement (TBJ- Unlike these dipole oscillations, for a finite field the atom en- 
circles the center of the parabolic lattice with a frequency defined by both the strength of 
the parabolic confinement and the magnetic field magnitude. This 'encircling frequency' is 
shown to be well seen in the quantum wave-packet dynamics. Thus the artificial magnetic 
field can be measured by measuring the encircling frequency. 

* * * 
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